Instability of Anisotropic Fermi Surfaces in Two Dimensions 
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The effect of strong anisotropy on the Fermi line of a system of correlated electrons is studied in 
two space dimensions, using renormalization group techniques. Inflection points change the scaling 
exponents of the couplings, enhancing the instabilities of the system. They increase the critical 
dimension for non Fermi liquid behavior, from 1 to 3/2. Assuming that, in the absence of nesting, 
the dominant instability is towards a superconducting ground state, simple rules to discern between 
d-wave and extended s-wave symmetry of the order parameter are given. 
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Since the appearance of high-T c materials the quest for non- Fermi liquid (NFL) behavior in two space dimensions 
has attracted great interest. The emergence of the quantum field theory approach to the study of electronic systems 
[1| seems to put severe constraints on the possibility of NFL behavior, at least for continuum models. Following Ref. 
[1| no NFL behavior is to be expected in D=2 provided that the coupling is weak, the interaction is short-ranged, 
and the Fermi surface is isotropic. On the other hand, Anderson [0 and coworkers support strongly the opinion that 
Fermi liquid behavior will not occur in D=2 under quite general assumptions, which make use of analogies to ID 
behavior. The basis of Anderson's approach is the discrete 2D lattice. On the other hand, it has been proved that 
the critical dimension for Luttinger liquid behavior in isotropic systems is D = 1 [£1. 

This paper addresses the role of anisotropy of the Fermi surface in inducing instabilities of the Fermi liquid fixed 
point in 2D. The cases of an extreme anisotropy due to either nesting or a Van Hove singularity are known to produce 
drastic changes in the system . We study here the effect that a strong anisotropy of the Fermi line can have in a 
continuum model in the absence of nesting and away from a Van Hove singularity. 

In between the case of a circular Fermi surface and the existence of a Van Hove singularity lies the possibility that 
the Fermi surface has inflection points which separate portions where the curvature has opposite signs, as shown in 
Fig. |j This situation takes place for a finite range of fillings, and requires no precise fine tuning of the number of 
electrons in the system. In particular, for the t — t' Hubbard model in a square lattice, inflection points at the Fermi 
surface are present for fillings below the Van Hove point ( — 8t' + 16i' 3 /t 2 < Sf < — 4t', where t < 0, t' > 0). 

Scattering processes which involve two opposite inflection points are analogous to the enhanced scattering present in 
nested Fermi surfaces Q. It can be shown that nesting, in any dimension, is a marginal perturbation. The correction 
that nesting induces to the bare interaction (in dimcnsionless units) is of order A log A, where A is the high-energy 
cutoff jjj . A Van Hove singularity induces similar effects in 2D || . 

In order to study the effects induced by inflection points at the Fermi surface, we first analyze the electron-hole 
polarizability (the bubble in Fig. |2](a)) for a wavevector Q which connects two such points at opposite ends of the 
center of the Brillouin Zone (BZ). The imaginary part of the diagram, Imn(Q, uj), measures the density of electron-hole 
pairs with energy oj. In the case of perfect nesting, each strip of holes of energy oj/2 below the Fermi level coincides 
with a strip of electrons at energy uj/2 above the Fermi level. Hence, Imn(Q,w) ~ M°, and the corresponding 
coupling is marginal. Near an inflection point, the states of equal energy no longer fill a straight strip in momentum 
space. Using a reference frame such that the x-axis is parallel to the Fermi surface, we can write: 

e(q) w v F q y + aq%, + (3q^ 
e(q + Q) w -v F q y - aql+ pq* (1) 

The number of electron-hole pairs of energy uj available are given by the region in momentum space which satisfy: 
u) — did < e e (q) — £/i(q + Q) < w. 

Setting q y ss lo/(2vf), these restrictions imply that the region available is bounded, in the transverse direction, by 
q-x oc q\J^ oc w 1 / 4 (note that the cubic term does not break the correspondence between the equal energy slices at q and 
q + Q). This implies that Imn(Q, lu) oc |oj | x / 4 . The previous argument can easily be extended to an arbitrary number 
of dimensions. An inflection point is characterized by a dispersion relation which is linear in the direction normal to 
the Fermi surface, cubic in one of the transverse directions, and quadratic in all the others. In D dimensions, this 
leads to: 
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ImII(Q,a;) oc \uj\ — +~ (2) 

Hence, the critical dimension is D c = 3/2. The same argument, when applied to the standard 2kp scattering which 
connects two opposite points of an isotropic Fermi surface, gives ImII(2/ci?,w) oc \uj\^ d ~ 1 ^ 2 , and the corresponding 
critical dimension is D c = 1, in agreement with Ref. Q. 

The preceding analysis shows that inflection points enhance the instabilities of the model. Among the possible 
broken symmetry ground states, we will assume that superconductivity prevails. The model has no tendency towards 
ferromagnetism, as the density of states needs not be specially large. A spin density wave, or a charge density wave, of 
momentum Q is possible. The associated gap, however, will only span a small fraction of the Fermi surface, centered 
around the inflection points. Moreover, only two of the eight possible points of a square lattice will take part in 
the symmetry breaking. As it is shown below, the system is unstable against anisotropic superconductivity. This 
instability opens a gap over large regions of the Fermi surface (except at the inflection points). Thus, this instability 
is a more robust alternative than a spin or a charge density wave. 

The possible processes which renormalize the particle-particle scattering are shown in Fig. ^. In a 3D isotropic 
system with repulsive interactions, they lead to the Kohn-Luttinger instability [To|| , which gives rise to superconduc- 
tivity at high angular momenta. Note that this instability is absent in 2D isotropic systems fiUM. For short-range 
interactions, the first and second diagrams in Fig. |2|cancel. This fact implies that pairing, if present, cannot arise from 
the propagator of a quasiparticle, in a similar fashion to phonon mediated pairing. In that case, superconductivity 
can be expressed in terms of a diagram like the first one in Fig. 0. If the electrons to be paired and the electrons 
which give rise to the quasiparticle are the same, diagrams like the first one in Fig. |^ are cancelled, to all orders, by 
diagrams like the following one in the same figure. 

In general, we need to consider the scattering of a pair of electrons of momenta k, — k into p, —p. We label these 
processes by the angle tfi between k and the x-axis, and p and the x-axis, <\>' . In the following, we will label each 
point at the Fermi surface by the angular variable, cj>. The value of the corresponding diagram in Fig. ^|(c) is denoted 

n(0,0')- 

We first analyze contributions of the type II(</>, <f). We assume a local electron-electron interaction, U . Then, using 
Eq. (|l|) we have 



llU,cb) =-i 1 



2tt) 3 



J dtu J dV~ e(q)/A e-< q+Q + 5Q >/ A C?(a;, q)G(w, q + Q + SQ) (3) 



where Q + <5Q is the vector which connects the point at the Fermi surface with momentum k with that at — k. The 
parameter A is a smooth cutoff in energies which is used to implement the renormalization group (RG) scheme. 

Using the new variables: e = e(q — <5Q/2) and e = e(q + Q + <5Q/2), the real part of H(4>, (f>) becomes, for the 
convex part of the Fermi line, 
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where the only dependence on comes through the curvarture: / = 12/3 SQ 2 — 6a SQ X . A similar expression holds 
for the concave part of the Fermi surface, with the same respective limit behaviors at large and small values of the 
cutoff A. 

The result (@) can be rescaled as a function of the dimensionless parameter - ^ , what allows us to study the 
relative intensity of the corrections in different zones of the Fermi line. We sec that, close to the inflection points 
(/ -C A), the vertex function reaches maximum values 

A l/4 

RellO^)--^ (5) 



while for points with / 3> A 



A 1 / 2 

Ren(0,0)^^ (6) 
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Equation (g) reproduces the 2/cf behavior of the response function of an isotropic Fermi liquid. We observe that 
the function 11(0, 0) has minima at points where the value of the curvature is a maximum. Its natural period is 
determined by the symmetry of the Fermi line, which in a square lattice is in general tt/2. 

By adding 11(0, 0') to the bare particle-particle scattering, U, we obtain the full irreducible vertex in the Cooper 
channel. This vertex, when inserted into ladder diagrams like the ones depicted in Fig. ||, gives the possible su- 
perconducting instabilities. The ladder summation implicit in Fig. || respects the lattice symmetries. We write the 
irreducible vertex as: 

U + ll(0i , 2 ) = c + ci (cos 0i cos 02 + sin 0i sin 2 ) 

+C2 cos 20i cos 202 + C3 sin 20i sin 20 2 

+C4cos40i cos 402 + C5 sin 40i sin 402 + ... (7) 



The effective interaction in the Cooper channel r(0i,02) has a similar expansion, with coefficients r„ whose cutoff 
dependence is given by: T n = c„/(l — c n /(47r) log(A/Ao)). Pairing instabilities are associated to negative values of 
the c„'s in Eq. (0). 

Let us now assume that the expansion of 11(0, 0') in angular harmonics can be truncated at fourth order, as shown 
in Eq. (?]). Then, from Eqs. @ and @, we can write: 11(0,0) m a + b sin 2 40, where: a ~ 0{K 1/2 / f 1/2 ) and 
b ~ 0{A 1 ' i /(3 1 ' A ), so that: 

co + ci + c 2 + c 4 - U ~ Oik 1 ' 2 /! 1 ' 2 ) 
c 3 - c 2 ~ 

c 5 - c 4 ~ 0(A 1 ' 4 /P 1 ' i ) (8) 

In order to close this system of equations, we compute II(— 0, 0). This function is ~ 0(A 1 ' 2 / f 1 ' 2 ) if = nn and 
shows a crossover to a behavior ~ O(A) for other values of 0. Hence, 11(0, —0) w c + d cos 2 0, where: c ~ 0(A) and 
d ~ 0(A 1 ' 2 / f 1 ' 2 ). It implies that: 

c - ci - c 3 - c 5 - U ~ 0(A) 

2ci ~ 0{A 1 ' 2 / f 1 ' 2 ) 

C2 + C 3 ~ 

c 4 + c 5 ~ (9) 

From (JsJ) and (|9|) we obtain that the dominant modes in the limit A — > are c - 17 ~ -c 4 ~ c 5 ~ ©(A 1 / 4 //? 174 )- 
Hence, C4 is the only negative coefficient to that leading order and the system has an instability towards extended 
s-wave pairing, with a gap which has two nodes in each quadrant of the BZ. We stress that this argument follows the 
standard BCS procedure, reformulated in RG language. The only approximation involved is the truncation of the 
expansion (Q) to fourth order. 

We can extend the previous analysis to anisotropic Fermi surfaces without inflection points, like the round shaped 
one shown in Fig. |], at the other side of the Van Hove filling. Then, the vertex n(0, 0) is always ~ 0(A 1 ' 2 / f 1 ' 2 ) but 
it is modulated by the local curvature of the Fermi surface. Hence, 11(0,0) has maxima for = 7r/4 and equivalent 
points, where the curvature reaches its minimum value. The same argument gives that the minima of n(0, 0) lie at 
= and equivalent points. The vertex is modulated by n(0, 0) w a + b sin 2 20, where: a ~ b ~ 0(A 1 ' 2 / f 1 ' 2 ). A 
similar analysis to that applied above leads to the conclusion that C2 is now the dominant coefficient that becomes 
negative, implying a superconducting instability with d x 2_ y i symmetry, in agreement with previous calculations for 
Fermi surfaces of similar shapes jL2j. 

We have assumed that the shape of the Fermi surface does not change much upon lowering the cutoff, A. This is 
the case if only the electrons near the Fermi surface are affected by the interaction, that is, A ~ U <C Ep, where Ep 
is the Fermi energy. Different superconducting instabilities may arise in dilute systems where this approximation is 
not valid . 

In conclusion, we have extended the standard renormalization group treatment of interacting electrons to systems 
with strongly anisotropic Fermi surfaces. We have shown that the existence of inflection points change the scaling 
properties of the interactions. The critical dimension which signals non-Fermi liquid behavior is 3/2. The supercon- 
ducting vertex, near the inflection points, also changes its scaling properties. This feature allows us to propose a 
simple scheme to analyze the possible superconducting instabilities of the system (the Kohn-Luttinger mechanism). 
We have extended the analysis to other anisotropic Fermi surfaces in the square lattice. The superconducting order 
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parameter is determined by the modulation of the curvature around the Fermi surface. When the curvature has two 
zeros in each quadrant of the BZ, we find extended s-wave pairing. The gap has two nodes in each quadrant of the BZ. 
If the curvature shows only one minimum, we obtain d-wave pairing, which can be x 2 — y 2 , if the region of minimum 
curvature is along the diagonal of the BZ, or of the xy type, if the minimum of curvature is along the axes. In all 
cases, the gap has a maximum in the region of maximum curvature of the Fermi surface. This region has the lowest 
density of states. It is interesting to note that a screening diagram, like the first one in Fig. |2|, favors a gap which 
spans regions of high density of states. An unusual feature of the Kohn-Luttinger instability is that it is induced by 
an antiscreening diagram, the third in Fig. ^|. 

For the particular case of the t — t' Hubbard model, we propose a phase diagram with extended s pairing for fillings 
below the Van Hove singularity, and d x 2_ y 2 for fillings above it. At the Van Hove filling the system shows magnetic 
instabilities. Different approximations indicate that the ferromagnetic instability is dominant for a finite range of 
parameters |l3[ ]. For values of t' below the ferromagnetic regime, though, a finite window for superconductivity with 
d x 2_ y 2 order parameter develops, as confirmed by a numerical diagonalization approach E3]. 

The closeness of a transition towards a d-wave and an extended s-wave superconductor found here implies that, in 
the absence of perfect tetragonal symmetry, a mixture of the two is likel y Il5[| . Such a possibility is consistent with 
recent photoemission experiments in the overdoped regime of BiSCCO |16| . Alternative scenarios for the interplay 
between the d-wave and extended s-wave symmetry of the order parameter have been proposed invoking different 
effective interactions [jl7| The framework presented in this Letter, on the other hand, links in a natural fashion 
the symmetry of the order parameter to the shape of the Fermi surface. Finally, our study can be relevant for the 
2D anisotropic superconductor Si'2Ru04 |20|, which is well into the weak coupling regime and has a superconducting 
state that influences a narrow energy range near the Fermi level (see, however, Ref. 
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FIG. 1. Different shapes of the Fermi line for the t — t' Hubbard model about the Van Hove filling. Two opposite inflection 
points are marked on the figure. 
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FIG. 2. Total interaction vertex in the BCS channel. 



FIG. 3. Ladder summation in the BCS channel used to describe the possible superconducting instabilities of the system. 
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